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The relaxation to equilibrium in many systems which show 
strange kinetics is described by fractional Fokker-Planck equa- 
tions (FFPEs). These can be considered as phenomenological 
equations of linear nonequilibrium theory. We show that the 
FFPEs describe the system whose noise in equilibrium funfills 
the Nyquist theorem. Moreover, we show that for subdiffusive 
dynamics the solutions of the corresponding FFPEs are prob- 
ability densities for all cases where the solutions of normal 
Fokker-Planck equation (with the same Fokker-Planck opera- 
tor and with the same initial and boundary conditions) exist. 
The solutions of the FFPEs for superdiffusive dynamics are 
not always probability densities. This fact means only that 
the corresponding kinetic coefficients are incompatible with 
each other and with the initial conditions. 
PACS No. 05.40.-a, 05.70.-a, 02.50.-r 



I. INTRODUCTION 

Different physical systems, like polymers chains, mem- 
branes, networks and other generalized Gaussian struc- 
tures, often show a long temporal memory due to the 
complex hierarchical organization of the modes of their 
motion. On the other hand, the response of these sys- 
tems to an external perturbation stays linear for a wide 
range of parameters p]-0]. As recently suggested, the 
response dynamics is well-described by dynamical equa- 
tions introducing fractional time-derivatives instead of 
whole- number ones From the thermodynamical 

point of view, the systems do not show any peculiarities 
close to equilibrium in contact with a classical heat bath. 

The relaxation to equilibrium in such systems is thus 
described by fractional Fokker-Planck equations (FF- 
PEs), which follow as phenomenological linear response 
equations. The corresponding equations are especially 
opular in application to a slow (subdiffusive) dynamics 
and where introduced ad hoc much before the micro- 
scopic basis for such equations got clear. 

We show that the typical FFPEs with the fractional 
derivative in front of the normal Fokker-Planck operator, 
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are the only possible variant for description of nearly 
equilibrium systems showing linear response, since they 
(and only they) fulfill the Nyquist theorem which con- 
nects linear response behavior with the noise spectrum 
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The value of 7 = 1 corresponds to an identity transfor- 
mation, leading to the case of pure diffusive behavior, 
the cases with 7 < f correspond to subdiffusive behav- 
ior, and the case 7 > 1 to a superdiffusive dynamics, like 
one considered in Ref. ||. 

Other forms of FFPEs are known, e.g. a Galilean- 
invariant form [j5]|8| which appears quite naturally when 
describing transport in a given velocity field, and forms 
with different fractional time-derivatives in front of the 
first and second spatial derivatives, which may appear as 
dynamical equations in many other contexts (economics, 
biology, etc.). They do not apply to cases of thermody- 
namical relaxation close to equilibrium. 

We show that FFPEs which describe subdiffusive dy- 
namics always have thermodynamically sound solutions 
when the corresponding normal Fokker-Planck equation 
also has them. Such solutions are subordinated to the so- 
lution of a normal Fokker-Planck equation with the same 
initial/boundary conditions. The situation with the su- 
perdiffusive dynamics is different: here not all combina- 
tions of external potential, diffusion coefficient and mem- 
ory kernel give rise to physical solutions (positive prob- 
ability densities), as it is e.g. the case for a fractional 
generalization of diffusion with drift. We discuss why 
it is so and exemplify this situation by processes subor- 
dinated to the solutions of a generic transport equation 
(related to a Liouville equation). 



II. FFPE'S AS A PHENOMENOLOGICAL 
LINEAR RESPONSE THEORY 

Let us first discuss the properties of FFPEs as phe- 
nomenological equations being very similar to nor- 
mal, diffusive Fokker-Planck equation (FPE). Within 
standard phenomenological linear nonequilibrium theory 
|)|.[f0[ the diffusion equation in a weak external field (i.e. 
a forward Fokker-Planck equation) follows as a conse- 
quence of local conservation of probability, 
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and a phenomenological linear response assumption 
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j = A (1) fP - \Wgrad P, 
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where A'- 1 -' = fi and A*- 2 -* = D are the kinetic coefficients 
(the mobility, and the diffusion coefficient, respectively). 
The phenomenological interpretation of the second equa- 
tion is that the current in our system can be caused by 
weak external field (and follows the Ohm's law) and by 
concentration gradient (the first Fick's law), and that 
both effects are independent as long as deviations from 
equilibrium are small. 

In general, the linear response can be retarded and 
then follows the equation 



j(t) = {f(t')P(t')} - $f } {grad P(t')} . 
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Here $t are typically causal integrals of convolution type: 
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where the lower integration limit to can be either finite or 
infinite. Here we again assume behaviors typical for the 
systems close to equilibrium. Inserting Eq. (|J) into Eq. (||) 
we get a nonmarkovian (nonlocal in time) Fokker-Planck 
equation of the from 
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(here we restrict ourselves to a one-dimensional case). 

Evaluating the first moment Mi (i) of the distribution 
P(x, t) under influence of a homogeneous force / we get 
that the evolution of the response follows the equation 
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from which it is clear that the operator is exactly 
the one describing the linear response of the system, so 
that the inverse operator corresponds to the system's 
impedance. 

Let us consider the noise produced by our system at 
equilibrium. The fact that the system is equilibrated 
means that it was created long ago, so that to — > — oo. 
Let us consider a Green's function of the equilibrium sys- 
tem, fulfilling the equation: 



d_ 

di 



G(x,t) = -& t 



d_ 

dx 

dx 2 



fG(x,t) 
G(x,t) 



(9) 



S(x)6(t). 



The Fourier-transform of the Green's function in both 
spatial and temporal domain is given by: 
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having a solution 
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Now, we are interested in the power spectrum of the equi- 
librium (/ = 0) noise generated by our system. Let us 
consider the second moment of G in frequency domain, 
M 2 (u) = G(k,u)\ k=0 = 2$( 2 Ww 2 . Note that x 
is the time-integral of the instantaneous velocity, so that 
the power spectrum of velocity (current) is exactly 
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Note also that the noise at equilibrium fulfils the Nyquist 
theorem [JO] , according to which 

S v (cj) = 2fc B TRe$ (1) (w), (13) 

so that the operators $W and $^ 2 ' are not independent: 

Re$W(w) = fc B TRe$( 2 )(w), (14) 

which for ^-operators of the fractional derivative type 
imply that $( 2 ) = ksTQ^. All equations with <f>( 2 ) = 
/cbT<I> ( 4) are thermodynamically sound: they fulfill the 
generalized Einstein relation and describe the relaxation 
to a Boltzmann distribution, which properties follow also 
from the microscopic description of the corresponding 
generalized Gaussian structures [Q. The equations with 
independent <&W and & 2 ^ will typically lead to behav- 
ior at variance with predictions of equilibrium thermo- 
dynamics. 

Note that the most systems for which the fractional 
dynamics was applied are "normal" although complex 
situations like polymers, membranes or fractal webs. In 
what follows we discuss only the case which describes 
such systems close to thermal equilibrium, for which the 
generalizations of FPE like Eq.([l]) can be considered as 
thermodynamically sound phenomenological laws. We 
also note that equations like Eq. ([!]) can be derived within 
the framework of stochastic approach |Q] , where they ap- 
ply to situations close to thermal equilibrium. On the 
other hand, the equations with different temporal oper- 
ators are also widely used: an example is a Galilean in- 
variant FFPE of Ref. ||. This equation appears quite 
naturally when describing transport in a given veloc- 
ity field, i.e. when our system is in a contact with a 
strongly nonequilibrium flow of fluid (a river instead of a 
bath!). Other variants with different orders of fractional 
temporal derivatives may appear as dynamical equations 
in many other contexts (economics, biology, etc.) but 
would never apply to the case of thermodynamical relax- 
ation in a system close to equilibrium, since they violate 
Eq. ( |l4[) . The situation with the systems whose dynamics 
shows linear response but is described by the FFPEs of a 
type different from one considered above is similar to one 
which arises when negative temperatures are considered 
pjfl: the systems described by such dynamics can live 
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as isolated systems but can not be in equilibrium with 
any " normal" macroscopic bath. Interacting with a heat 
bath, such systems will gain or lose energy until they 
leave the linear response regime and get a noise spec- 
trum conformal with equilibrium (and with a Boltzmann 
energy distribution). 

Note that the FFPEs like Eq. (jl|) were proposed for 
processes with finite increments (like continuous time 
random walk processes) or ones with continuous trajec- 
tories (fractional Brownian motion), situations for which 
the assumption of the local (differential) conservation law 
is proved. The related thermodynamical considerations 
show that a system whose noise does not posses any sec- 
ond moment (Levy- noise) , does not fulfill local conser- 
vation, Eq.(^), and can hardly exhibit linear response, 
a fact found in Ref. Jl3| on an example of subordinated 
processes. This case is addressed in Ref. ]l4| and leads 
to a different form of FFPEs with fractional spatial op- 
erators. 

Considerations based on linear noncquilibrium thermo- 
dynamics are somewhat too general, since Eqs.(g) and 
(||) guarantee the overall conservation of the value P but 
not the fact that this P is a nonnegative quantity. The 
same equations will apply for electric charge and cur- 
rent (which can be of both signs and may oscillate) and 
for density or temperature, which are essentially nonneg- 
ative. Thus, in order to check that the corresponding 
equation is thermodynamically sound one has to prove 
that if the initial condition corresponds to a nonnegative 
density P(a:,0), the density P(x,t) will stay nonnegative 
during all the following evolution. Since we concentrate 
here on the properties of relaxation to equilibrium, the 
force term and the diffusion coefficient in our system will 
be considered time-independent. 

The proof of the non-negativity of solution for the 
force-free case was given in jL5| for the subdiffusive case. 
We show that the same is the case for the arbitrary exter- 
nal force. Namely, we shall show that all solutions of FF- 
PEs with 7 < 1 in arbitrary time-independent potential 
force field are thermodynamically sound, and describe 
the transport of a positively defined density. Moreover, 
we show that superdiffusive equations with 1 < 7 < 2 do 
not always possesses physically sound solutions, unless 
some additional conditions are fulfilled. Fokker-Planck 
equations of the type of Eq.(|l]) with 7 > 2 seem to con- 
tradict physical sense. However, the superballistic be- 
havior (say, Levy flights) can be described by the FFPEs 
of a different class, see Ref. 



III. THE SUBDIFFUSIVE CASE: TEMPORAL 
SUBORDINATION 

Let us first consider the subdiffusion case, < 7 < 1. 
Note that the solution of subdiffusive FFPE under time- 
independent force can be put in the following form: 
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with yC^/r 1 ' 7 , 7, —7) being an extreme (one-sided) Levy- 
stable law of index 7 |5p6[], and F(x, r) is the solution of 
" normal" FPE under the same force and the same initial 
conditions: 

^F(x, t) = -pJLf{ x )F{x, t) + k B T^F(x, t). (17) 

To check this let us take the Laplace-transform of both 
sides of Eq.(|l|), and note that this transform acts only on 
the t variable, which appears in the Eq. (|l5[) as a param- 
eter: The Laplace transform in t of Eq.dlS) reads: 
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The fractional temporal differentiation leads to 
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Note that for t — > the T-functions are strongly concen- 
trated, so that T(T,t) -> S(t) and P(x,+0) = F(x,+0). 
Changing now to a new variable A = u 1 we recognize in 
Eq.(|l9|) the Laplace-transform of the "normal" FPE with 
the same time-independent force and the same initial con- 
ditions. Thus, the solution of FFPE can be obtained from 
the solution of FPE by immediate integration. Moreover, 
each functional of such solution (e.g. any moment) can 
be immediately obtained by weighing the corresponding 
functional of the FPE solution with a probability dis- 
tribution, Eq. (|l6|) for which useful analytic representa- 
tions are known. Thus, the equations with 7 < 1 in any 
(temporally constant) force field / obey regular Boltz- 
mann thermodynamics and correspond to the transport 
of a positively defined density. Our result generalizes 
the mathematical treatment of Schneider and Wyss and 
shows that the solution of a FFPE describing subdiffu- 
sive transport in external potential is a probability den- 
sity whenever the solution of a normal FPE in the same 
potential is one. The generalization to higher dimensions 
is evident. Note that our discussion here parallels that of 
Ref. |l(J where the fractional Kramers equation is con- 
sidered. 

Note that Eq. ( |l9| ) shows an extremely interesting prop- 
erty of free relaxation of the systems described by subdif- 
fusive FFPEs, namely the fact that the solution of Eq.(Q) 
having a form of convolution (linear response with a long- 
time memory kernel) can be represented in the form of 
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subordination, i.e. they correspond to the behavior of 
the system whose development is governed by its own 
internal clock, which is not synchronized with our phys- 
ical time |l9|| . The first reasonable use of this fact can 
be probably attributed to P.K. Clark, see Ref. |2(| for 
the discussion of the role of subordination in economical 
processes. This (operational) time is a variable which 
is Laplace conjugated to u 7 , and can be considered as a 
real time variable, since it is monotonously growing in our 
physical time and allows to order the events sequentially. 



IV. SOME PROPERTIES OF TEMPORAL 
SUBORDINATION 

The integral transform, Eq.([l5|) will be called a subor- 
dination transformation (ST), the term "time-expanding 
transformation" (TET) will be reserved for those with 
7 < 1. In order to sharpen the instruments needed for 
understanding the consequences of Eq. ([!]) let us discuss 
some properties of STs with T(r, i)-functions from the 
class discussed above. The physical time t will be called 
the outer variable of the function T(r,t), and the vari- 
able r (operational time) over which the integration is 
performed will be called its inner variable. Note that the 
STs are just a type of transformation typically arising in 
a context of separation of variables (the eigenfunctions 
decomposition with integration or summation over the 
eigenfunctions numbered by the eigenvalue r). In the 
eigenfunction decomposition we start from the solution 
of a time-independent equation (zero order in time) and 
rise the order of temporal derivative (say to one) by ap- 
plying an equation of the type of Eq.(15). In the case of 
subdiffusive FFPE we proceed in the opposite direction: 
we lower the order of a temporal derivative by applying 
subordination. As we proceed to show, a superposition 
of two functions of the type of T(t, t) (with indices 71 
and 72) discussed above, is again a function of the same 
class with the index 7 = 7172. Let us suppose that both 
indices, 71 and 72 are less than unity. We know that the 
Laplace-transform of T(r, t) in its outer variable reads: 



T(r, u) 



,7-1 



exp( 



We thus get: 



inner variable: they are nonnegative and integrable. On 
the other hand, the T-functions rising the order of the 
temporal variable have a Laplace-transform in the outer 
variable which reads: 



(20) 

which is again a Laplace-transform of a T-function in its 
outer variable. Using this fact once again we get 

/>oo 

T*{t,u) = u^ 1 / T^OeM-Z^H 
Jo 

= u 7271 " 1 exp(-™ 7172 ). (21) 

Thus, parallel to the Levy-case of Ref. fl3[| , the super- 
position of two TETs is a TET again. Note that all 
T-functions with 7 < 1 are probability densities in their 
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1 exp(-Tu 1 / 7 ), (22) 



i.e. belong to the same class of functions than T's 
themselves, but with 7* = 1/7 > 1. Note that the 
transforms T 7 and T~ = Tl/ 7 , lowering and rising the 
order of the FFPE to the same amount are the in- 
verse of each other: the Laplace-transform of T^T~ is 
cxp(— ut) so that T 7 T~ corresponds to a 6(r — t). More- 
over, we have to note that the r-integral of T(r, t), 
N(t) = J T(r, i)dr being an inverse Laplace-transform 
of N(u) — J °° u 7_1 exp(— Tu f )dr = u~ 1 , is equal to 1 
both for TETs (7 < 1) and inverse (7 > 1) transforms, 
so that both the subordination and the inverse trans- 
formation keep the overall normalization of the possible 
PDFss as functions of coordinates. 

The T~-functions are not PDFss of r since they may 
take negative or even complex values. Let us fix some 
value of t and consider the limiting value of the integral 
I(t) = f™T(r,t)dt, which can be expressed in terms 
of T(t,u): I(t) = \im u ^ (u1- 1 e~ TU ' , ). For 7 < 1 the 
corresponding integral diverges being positive. On the 
other hand, for 7 > 1 I(t) = 0, which means that 
the function T(r, t) either changes its sign or vanishes 
identically. The last is not the case since the integral 

still di- 



7-1, 



h{r)= f °°tT(r,t)dt = -±(u 

verges for 1 < 7 < 2 (for larger values of 7 the integrals 
J„( T ) = f Q °° t n T(r, t)dt with n > 7 - 1 still diverge). 



V. THE SUPERDIFFUSIVE CASE 

Our derivation of the FFPE and its formal solution 
through subordination are valid independently on the 
particular value of 7. The fact that T 7 (t, t) for 7 > 1 is 
not nonnegative does not mean that the integral Eq. ( |T5| ) 
takes negative values: it solely means that the non- 
negativity of the physical solutions of FFPEs does not 
follow from the non-negativity of the physical solutions 
of the Fokker-Planck equation, and that the variable r 
can be no more interpreted as an internal time governing 
the system's evolution. On the other hand, Eqs.([l^) and 
( |l9| ) are still valid as a representation of a formal solution 
of the FFPE. We shall refer to such formal solution as 
following from a pseudo-subordination. In some special 
cases of pseudo-subordination one can still can guarantee 
that the corresponding solution is a probability distribu- 
tion, as it is e.g. the case for force-free transport for 
7 < 2, in other cases the solutions are not PDFs as it is 
e.g. for 7 > 2. 



4 



A. Pure superdiffusion: Relation to a wave equation 

Let us consider a purely diffusive situation without ex- 
ternal force, 



(23) 



The Laplace-transform of this function in i-variable 
reads: 



F(x,u) = -u 1/2 exp (- \x\ \/u) 



(24) 



Let us now use Eq.(|l^), and get P(x,u) for arbitrary 7: 

P(x, u) = v?- 1 F{x, v?) = \u> /2 ~ 1 exp (- \x\ u 7/2 ) . 

(25) 

The function P(x, u) belongs to the class of functions 
T(r, t) given by Eq. (|l6|) , but with change of 7 to 7/2: 



P(x,t) = -T l/2 (\x\,t). 



(26) 



Note that Eq.(|26|) gives the representation of the su- 
perdiffusive propagators in terms of the Levy-functions, 
which simplifies the general result of Ref. |J. Since we 
know that T 7 (£, t) is a positive function of its both vari- 
ables for t > and 7 < 1, in the case of free propagation 
P(x,t) is positive for all 7 < 2. The case 7 = 2 corre- 
sponding to 



P(x,t) = ±6(\x\-t) 



(27) 



describes a special case of ballistic propagation. The 
equations of index 7 > 2 (describing a process which 
is faster than ballistic one) do not correspond to trans- 
port of positive probabilities, since the functions T 7 with 
7 > 2 are no more non-negative. 

We have seen that although the non-negativity of 
the solution is not mathematically guaranteed by the 
FFPE with 7 > 1 itself, the equation still can pos- 
sess physically reasonable positive solutions describing 
superdiffusive transport. Let us discuss now the reasons, 
why it is so. Let us note that the solution P(x, it) = 
it" 1 / 2 exp (— |x| y/u) itself can be considered as subordi- 
nated to a process described by Eq.(p7|) (corresponding 
to ^'(x, u) — \pK exp (— \x\ u) ) with a " subdiffusive" sub- 
ordination function Tx/ 2 (t, t), so that the whole process 
can be considered as a superposition of two subordina- 
tion transformations, leading to the overall behavior with 
7* = 7/2. The process subordinated to a ^-functional 
form under operational time given by T- 7 / 2 (r, t) is, of 
course, exactly the solution given by Eq.(E9) discussed 
before. 

We note here that the two 5-pulses described by 
Eq.(p7|) are a solution of a wave equation (WE), 



d 2 * _ <9 2 * 



(28) 



The solution Eq. (|27|) is not a Green's function solution 
of a wave equation (known to be G(x,t) = \0{t — |x|) 
in one dimension, see Ref. p3[ ) but a solution cor- 
responding to a different initial condition, namely to 
G(x,t) — > 6(x)5'(t). The reason for this is easy to un- 
derstand: The limiting equation for the Green's function 
of a FFPE with 7 — > 2 is not a wave-equation, but a 
first-order integro-differential form, 



dF 
dt 



t pa 



d 2 F 
q dx 2 



dt' + 6{x)S(t), 



(29) 



which is obtained from a wave equation by temporal in- 
tegration. 



B. A problem of superdiffusion with drift 

Let us now consider processes being pseudo- 
subordinated to diffusive motion under time- independent 
homogeneous external force (i.e. the solutions of FFPEs 



dP 
~dt 



(30) 



with a = 1 — 7 < 0). The Laplace-transform of the 
corresponding Green's function of FPE, 



F(x,t) = 



1 



2\^Dt 



exp 



(x - fifty 
ADt 
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reads: 



F(x,u) = 
exp(^) 
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-v/(C 2 + «)A 2 



dt 
(32) 



where the variables A = xj2\fD and £ = fif ' /2yD (£ > 
0) are introduced (see 2.3.16.2 of Ref. p^]). Applying the 
T2-transformation to Eq.(32) we get: 



P 2 {x,u) = 



exp(2(A) 



2yfD y/C 



: exp 



-V(C 2 +« 2 )A 2 



(33) 



Let us show that Pz(x, u) is not a Laplace-transform of a 
probability distribution. Note that a Laplace-transform 
f(u) of any nonnegative function f(t) must be an abso- 
lutely monotone function, i.e. (— l) n -£-^f(u) > must 
hold for all u and n. To prove this it is easy to see that 
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— f(u) = I f(t)e- ut dt = I t n f(t)e- ut dt 



(34) 



On the other hand, the first it-derivative of P2{x,u), 
changes its sign at u being a root of the equation 



C 2 vV + C 2 )A 2 -«V + C 2 ) = 0. 



(35) 



The existence of positive roots of this equation for any 
(, =/= is clear since for u small the overall expression 
(whose sign is the same as the sign of -4-P^x, uj) is pos- 
itive, and for large u it is negative. Note that the func- 
tion gets to be absolutely monotone only when £ = 0, 
i.e. only in the case of free diffusion. This observation 
is of extreme importance since it shows that while the 
TETs (T 7 with 7 < 1), lowering the order of FFPE, al- 
ways lead to reasonable physical solutions, the inverse 
transformations, rising the order of FFPE do not always 
do so. 

Note that all functions P-)(x, t) obtained from diffusion 
with drift under pseudo-subordination are not probabil- 
ity distributions for all 7 > 1 . The Laplace-transform of 
the corresponding functions read 



P 1 {x, u) 



exp(2£A) u 

2Vd VC 2 



7-1 



■. exp 



V(C 2 + ^)A 2 



(36) 



The first derivative of P 1 (2, u) changes its sign at u being 
a positive root of 



((7-2K + (2 7 -2)C 2 )vV 
-7^(u 7 + C 2 ) = 0, 



C 2 )A 2 



(37) 



which function is positive for small u and negative for 
larger ones. Thus, the solutions of the FFPE of the type 
of a Fokker-Planck equation with 7 > 1 under homo- 
geneous, constant force are not probability distributions. 
Similar conclusions were drawn when considering the par- 
ticle's motion in a harmonic potential [E4J . 



VI. SUPERDIFFUSIVE CASE: SUBORDINATION 
TO A GENERIC TRANSPORT EQUATION 

In Sec. 5. 1 we have seen that the solutions of a diffusion 
equation (representing a behavior of a stochastic process) 
are subordinated to deterministic dynamics, described 
by a simple propagation of pulses with constant velocity 
and given by a wave equation. Is the wave equation (i.e. 
limiting superdiffusive FFPE with 7 = 2 for a force-free 
situation / = 0) a very special case, or there are some 
forms with / 7^ which still lead to reasonable solutions? 

It is clear that the second-order partial differential 
equation to whose solutions the solution of FFPEs could 
be subordinated would read: 



dx 2 



[B{x)P] 



(38) 



Eq.(38) includes the wave equation as a special case. 
Eq.(38) will be called the generic transport equation 
(GTE), and, parallel to a wave equation, has a dynamical 
(deterministic) nature. This equation (being a close rel- 
ative of Liouville equation) was considered by the author 
in a different context in Ref . |25[ : The GTE appears when 
restoring temporal dependence in a Pope-Ching equation 
for stationary random processes, Ref. |2qj . The meaning 
of prefactors here is: A{x) — (x(x)) and B{x) — (i 2 (x)), 
so that for a physical particle they are proportional to the 
acting force and to the particle's mean kinetic energy. 

Let us remind the procedure of derivation of GTE 
given in Ref. ^5j. The PDF of x, p x (x), is obtained 
as an ensemble-average (e.g. over the initial conditions) 
of the realizations for each of which 



p{x,t) = S(X(t) - x) 



(39) 



where X(t) represents the law of motion. The coarse- 
grained probability is then given by p(x) = (p(x, t)). Dif- 
ferentiating Eq.(B9) with respect to time one gets 



dp = _j c ^P_ = ^_ j V;j 

dt dx dx 



(40) 



since X is independent on x. Note that Eq.(f40|) is a 
Liouville equation, and the derivation here is parallel to 
one given in Ref. j2jj |. Applying the same procedure for 
the second time we get: 



d 2 p d d \ d 
dt 2 dx dt V / dx 



dx dt 
d_ 

dx 



dx 



Xp 



x 2 P 



(41) 



which equation is, of course, the as exact as the Liouville 
one. The GTE follows after the ensemble averaging, un- 
der which the corresponding conditional means appear 
instead of the instantaneous velocity and acceleration, so 
that Eq.@ reduces to Eq.@. Note that the GTE can 
be useless but is never false: its solutions describe all 
possible motions and are both dynamically and thermo- 
dynamically sound. These solutions are probability den- 
sities. On the other hand, the prefactors A and B arise as 
(nonequilibrium) ensemble averages and depend on what 
ensemble is used and thus on the initial conditions: a sim- 
ple example of this fact is considered below. The absence 
of the physical solution of Eq. (pq) means that the corre- 
sponding thermodynamical forces and kinetic coefficients 
defining A(x) and B{x) are incompatible with each other 
or with the initial conditions and would never appear as 
thermodynamical means. Moreover, even if the system 
as a whole is homogeneous and its physical properties 
are time-independent, the coefficients A(x, t) and B{x, t) 
can be time-dependent and will relax to the equilibrium 
values not faster than the distribution itself relaxes to 



G 



its equilibrium form, which explains the unphysical sort- 
time behavior of the solutions of superdiffusive FFPE in 
harmonic potential found in Ref. f24j . 

As an example of a processes subordinated to a solu- 
tion of GTE let us consider a simple oscillatory process 
taking place in the operational time of the system. The 
dynamic equation of the oscillator is 



-uix. 



(42) 



Let us consider the situation when the oscillator starts 
with zero velocity at x = —a so that A(x) = —lu 2 x, 
and B(x) = cu 2 (a 2 — x 2 ). Our process is described in 
operational time by a GTE 



d 2 F d 



dr 2 dx 



(^F) + ^[^(a 2 -x 2 )F] 



(43) 



with the initial conditions F(x,0) — 5(x + a), and 



OF{x,t) 



dr 



0, whose solution, as anticipated, reads 
f o cos cut) . Note that the coefficient B 



F(x, t) = °5(x 
depends explicitly on a, so that the form of equation de- 
pends on the initial energy of the oscillator. Equation 
( p3| ) is incompatible with any combination of initial con- 
ditions not leading to the same amplitude of oscillations, 
i.e. whenever one supposes lu 2 x 2 (0) + x 2 (0) ^ a 2 , and 
would lead in this case to negative or complex PDFs. 
The solution subordinated to F(x, r) reads: 



P(x, u) 



5(x + a cos lot)u' 1 exp(— Tu J )d,T 
1 



iiy/l — x 2 /a- 



: exp 



(~) 



exp( nu 1 ) 



(44) 



71=0 

■lf<- 



exp [i arccos (—-) u 1 ] 



Loa^/l-x 2 /a 2 l-exp(-^uT) 

For t — > oo (full dephasing) the corresponding PDF tends 
to a PDF to find an oscillating point at coordinate x. 
Thus, for u — ► 



1 



P(x,u-> 0) ~ - 



1 



u na^l - x 2 /a 2 



which corresponds to 
P{x, t — ► oo 



~k<x\J\ — x 2 /a 2 



(45) 



(46) 



a well-known solution for the invariant PDF for a clas- 
sical harmonic oscillator. Thus, the equations subordi- 
nated to GTE may describe partly coherent phenomena: 
their physical relation to a wave equation gets evident. 
It is interesting to mention that the Fokker-Planck equa- 
tion with the same coefficients as Eq.(ff3|) describes a har- 
monic oscillator in which the diffusion takes place in an 



inhomogeneous temperature field, T(x) ~ v/T- x 2 /a 2 , 
and such a diffusive solution stays physically sound both 
under subordination and under pseudo-subordination up 
to 7 = 2. 

The discussion above shows that for all A(x) and B(x) 
and initial conditions for which Eq.([38|) has real, non- 
negative, normalizable solutions, the free relaxation of a 
complex system under FFPE dynamics can be described 
as its deterministic development in its own (operational) 
time. The inverse (as we have proved on an example of a 
diffusion with drift) is not the case. Our considerations 
leave open the question whether all physically sound so- 
lutions of superdiffusive FFPEs are subordinated to ones 
of GTE, or this class is wider and includes some functions 
which are probability densities for 1 < 7 < 7* and cease 
to be probability densities for 7* < 7 < 2. 

It is important to stress that the fact that the solutions 
of superdiffusive FFPEs in general are not probability 
densities, does not devaluate the FFPEs as an instru- 
ment for description of complex relaxation phenomena, 
but shows that many combinations of thermodynami- 
cal forces, kinetic coefficients and memory functions will 
never appear as thermodynamic ensemble means. This 
means that the forces and kinetic coefficients can not be 
invented ad hoc, but must follow either from experiments 
or from microscopic considerations. In the case when the 
correct thermodynamical forces and the impedance of the 
system are known, its FFPE is uniquely determined. 

The consideration of GTE explains also our finding 
that the solutions for the force-free transport with 7 > 
2 are not non-negative. Such equations would describe 
processes subordinated to the solutions of the exact third- 
order transport equation. The exact equation with 7 = 3 
arising from applying a Liouville operator + ±Jj to 
the PDF P three times is a trinomial construct, with 
correlated coefficients in front of the first, the second, and 
the third spatial derivatives. This third-order equation 
is a generic form for transport equations of higher order. 
The equations subordinated to this one will have third- 
order structure in spatial variables and will be hardly a 
helpful tool, since they do not have any known classical 
counterpart whose solutions may be used for building 
new ones. 



VII. CONCLUSIONS 

Fractional Fokker-Planck equations (FFPE) with ad- 
ditional fractional time derivative in front of a normal 
Fokker-Planck operator appear within a usual linear- 
response theory when describing systems showing strange 
kinetics. We show that such form of FFPEs describes 
systems in a contact with a heat bath, since the noise 
in such systems in equilibrium (for t — > 00) fulfills the 
Nyquist theorem. Many other forms (e.g. with tem- 
poral fractional derivatives of different orders in front 
of first- and second-order spatial derivatives) are thus 
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ruled out as appropriate for describing situations close to 
equilibrium, although they may be appropriate for many 
other transport processes, as e.g. dispersion by flows. 
Using the fact that the solutions of subdiffusivc FFPEs 
with time-independent coefficients are subordinated to 
those of normal Fokker-Planck equations, we show that 
the FFPE solutions are probability densities in all cases 
when the usual Fokker-Planck equation has physical solu- 
tions. Thus, the free relaxation of a complex system un- 
der FFPE dynamics can be described as its development 
in its own (operational) time. The superdiffusive FFPEs 
do not posses physical solutions for arbitrary choice of 
force and diffusion coefficient. This does not devaluate 
the FFPEs as a tool for description of superdiffusive pro- 
cesses, but stresses the fact that the corresponding com- 
bination of thermodynamical forces and memory func- 
tions can never emerge as a thermodynamical ensemble 
average: the corresponding phenomenological equations 
have to be handled with care. 
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